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We examine the implications of a beat frequency modulated model of T Tauri accre- 
tion. In particular we show that measurements of the variability of accretion generated 
lines can be used in conjunction with existing photometry to obtain a measurement of 
the underlying photospheric and disc flux. This provides an independent way of check- 
ing spectral energy distribution modelling. In addition, we show how spectroscopy of T 
Tauri stars can reveal the inclination angle between the magnetic axis and the plane of 
the disc. 
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INTRODUCTION 



rffrp Classical T Tauri stars (CTTS) are believed to be ac- 
(ttfre ly accreting material from circumstellar disks, giving rise 
tfTp. pronounced veiling continuum, and characteristic emis- 
features such as H a and Ca II H + K. Periodic pho- 
tSlJietric variability is also often observed in these systems. 
' Jftefi standard picture is that this variability is caused by 
rqtational modulation of localised accretion shocks at the 
sQlar surface. This model cannot be applied in all cases, 
ljojyever, as many periods are inconsistent with the stellar 
rtreation as revealed by v sin i measurements (Bouvier et al 
lyS6; Bouvier et al 1993). Furthermore, there are a number 
of CTTS for which the period is seen to change (Bouvier et 
al 1994), an impossible situation if the period reflects the 
stellar rotation. In a recent paper (Smith et al 1994), we 
have discussed the possibility that the photometric period 
is caused by an interaction between a misaligned magneto- 
sphere and a clumpy disk, a model analogous to that de- 
veloped for horizontal branch quasi-periodic oscillations in 
low-mass X-ray binaries (Alpar & Shaham 1985). The pe- 
riod seen in the photometric monitoring would be the beat 
period between the stellar rotation period and the orbital 
period at the inner edge of the disk, or some harmonic of 
this period. Support for this model has been found in line 
profile variability of SU Aur (Johns & Basri 1995) which dis- 
plays evidence for periodic accretion with the same period 
as that observed in the photometry. 

In this paper we point out that if the periodicity is 
caused by beat frequency modulation, the amplitude of the 
variation depends on the non-axisymmetry of the magneto- 
sphere at the inner edge of the disc, which in turn depends 
solely on the inclination of the disk to the equatorial plane 



of the star. It is thus possible to recover the inclination of 
the disk plane from the photometry. In sections 2 and 3 we 
outline the model, in section 4 we show how the model can 
be applied in a specific case, and in section 5 we summarize 
our conclusions. 



2 THE MODEL 

The star accretes material from the inner edge of its disk as 
magnetic torques brake the disk material, which then falls 
in along the field lines. The rate of accretion will depend on 
the magnetic field strength, the surface density of the disk, 
and the exact nature of the interaction between the field and 
the disk. 

In the beat frequency model, the interaction between 
the non-axisymmetric field and inhomogeneous disk leads 
to variations in the accretion rate, with a characteristic fre- 
quency \io,nner — the beat frequency between the disc 
and the star (here io t „„ e r is the Keplerian frequency at the 
inner edge of the disk, and lu* is the spin frequency of the 
star). Assuming the clumpy structure of the disk remains 
constant for a large number of orbits (a necessary assump- 
tion if the T Tauri periodicities are to be explained in this 
way), the amplitude of the variations in the accretion lumi- 
nosity depends only on the difference between the minimum 
and maximum field strengths at the inner edge of the accre- 
tion disc. 

The rate at which material is stripped from the disk as 
a function of field strength is not well known. The torque 
exerted on a given disk annulus by a field which fully pene- 
trates the disk is proportional to the product of the poloidal 
and toroidal field components, B z and B<f, (Armitage 1995). 
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Here the poloidal field is assumed to be a pure dipole field, 
while the toroidal field is formed by distortion of the dipole 
field by the fast moving disk material. The magnitude of 
B,f, depends on which process is invoked to limit its growth. 
Here we follow Armitage (1995) and assume B<f, is limited 
by reconnection of field lines. In this case the equilibrium 
value of B,f, is of order B z . A realistic estimate of the torque 
is then 



Disc inclination as function ot Variability amplitude. 



T oc B^B Z oc B z 



(1) 



In the interests of generality, we use T oc B" z for our torque. 
We assume that the stripping rate and the accretion lumi- 
nosity are linearly dependent on the torque. Thus, the ratio 
of accretion luminosity at maximum field strength, f(Bmax) 
to that at minimum field strength, f(B mtn ), is 

f(B max ) (B 

max 
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Z 



(2) 



The measurement of Z can be used to gain a valuable insight 
into the geometry of a given T Tauri system, since it depends 
only on the ratio of maximum to minimum field strength at 
the inner edge of the disc. For a purely dipolar field, this 
ratio is; 



B„ 



B„ 



(1 + 3 sin 2 (9)" 



(3) 



(Lipunov 1987) We argued above that 7 = 2 for a field 
whose toroidal component is limited by reconnection. Thus, 
knowing Z, it is possible to recover the inclination angle 
between magnetic dipole axis and the plane of the disc, 9. 



3 PERIODIC ACCRETION AND EMISSION 
LINES 

By making the assumption that the dependance of accretion- 
generated lines on the field strength is the same as that of 
the continuum, we can use the ratio of maximum to mini- 
mum line flux in a number of lines to measure the ratio Z. 
Kuhi (1991) showed that, for a number of stars, a strong 
correlation exists between the flux in H a and the ultraviolet 
continuum, which is dominated by accretion luminosity, and 
also between the flux in Ca II K and the ultraviolet contin- 
uum. This provides support for our assumption above. The 
lines used should be high energy, so they are more likely to 
be driven purely by the violent accretion events at the star's 
surface, and should vary with the same period as the con- 
tinuum, which would demonstrate their link to the varying 
accretion. 

The periodic variability in the bolometric photometry 



AM = -2.51og 10 
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F + f{B m 



(4) 



Here AM is the amplitude of the photometric variability and 
Fo is the underlying, constant, combined luminosity from 
the star and disc. By making measurements of the maximum 
and minimum magnitude attained by the star in its periodic 
variations, and knowing Z from line measurements, we can 
solve for Fo and f(Bmtn), and hence find f(Bmax)- The 
equations for photometry in a band, X, are 
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Figure 1. The disc-magnetic field inclination as a function of 
the observed bolometric variability. Three different cases are con- 
sidered, with fmin equal to the star and disc luminosity, lOi^o, 
Fo, and O.li^o, corresponding roughly to accretion rates of 10 — 6 , 
10 and 10~ 8 Mat/r- 1 In all cases 7 is taken to be 2. 



and 



X + af(B„ 



(6) 



Where X m ax and X m tn are the maximum and minimum X 
band absolute fluxes, Xo is the X flux from the star plus 
the disc, and a is the fraction of the accretion luminosity 
appearing in the X band. Combining these equations gives 
us 
X 

m a. x -Xo 



z 



(7) 



X„ 
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(5) 



Xmtn — Xo 

Knowing these fluxes from the photometry, and Z from the 
line variability it is possible to solve for Xo. It is then pos- 
sible to solve for af(Bmtn)- 

Broadband photometric observations of periodic vari- 
ability have been published by a large number of authors 
(Bouvier et al 1993; Bouvier et al 1994; Attridge & Herbst 
1992). The amplitude of variability increases sharply in the 
B and U bands, possibly because the non-variable contribu- 
tion of the stellar photosphere is much lower here. Typical 
observed amplitudes are 0.5 magnitudes in U, decreasing to 
0.1 magnitudes in V, and less in R and I. 

Figure 1 shows the amplitude of the observed bolomet- 
ric variability produced by different field inclinations. The 
three curves represent different accretion rates, and hence 
different accretion luminosities. It can be seen that the am- 
plitudes produced are consistent with the observed photom- 
etry for relatively modest inclination angles. 

In addition to directly affecting the accretion and the 
emission lines involved, the inclined rotating magnetic field 
will also affect the surrounding disc and its emission. 

First of all, magnetospheric accretion removes the inner 
part of the disc up to the magnetospheric radius. With an in- 
clined dipole, the disc sees a varying magnetic field strength. 
Crudely translated into a varying magnetospheric radius, a 
dipolar field inclined at 30 degrees implies a variation in 
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Band Amp ^ F (%) DG Tau(MO) DG Tau(K7) 



Line A (A) EW min EW max Z 



Ca II H 


3968 


-27 


-60 


2.22 


Ca II 


8662 


-50 


-68 


1.36 


H a 


6561 


-75 


-105 


1.4 


Hb 


4861 


-24 


-42 


1.75 


He I 


5876 


-2 


-3.8 


1.9 



Table 1 . Equivalent width variations for five lines in the spectrum 
of DG Tau, from Guenther (1994). The implied vales of Z are 
given, assuming a one to one correlation between the equivalent 
width and line flux. 

the magnetospheric radius of 17%. Although this in practice 
won't translate directly into a change of the inner hole of the 
disc, as the disc will not adapt immediately to the change in 
the magnetic field, it will have effects on the emission from 
the inner part of the disc. Furthermore, the dissipation of 
the magnetic field in the disc will increase the emission in 
the near infra-red (Armitage & Clarke 1995; Kenyon, Yi & 
Hartmann 1995). Any asymmetries in the disc will then also 
produce variability in the near infra-red with the same pe- 
riod as the accretion. This might explain the near infra-red 
variability found in T Tauri stars on timescales from a few 
days to weeks (Hillenbrand, personal communication). 



4 AN ILLUSTRATIVE EXAMPLE 

Two sets of observations are necessary to complete the mea- 
surements we propose above. Firstly, a spectroscopic mon- 
itoring of high excitation energy accretion linked emission 
lines to determine Z, and secondly, a measurement of the 
photometric variability. 

Any emission lines observed would need to be linked to 
the accretion process. For this reason they should be seen to 
vary with the same period as the photometry. Such variabil- 
ity was observed by Guenther (1994), for 5 emission lines 
in three systems. Selecting high excitation energy lines is 
necessary, since such lines are expected to form in the hot 
accretion shocks, and are thus more likely to be purely ac- 
cretion driven. 

DG Tau is a CTTS believed to have an MO photospheric 
spectrum. It was one of the three stars studied by Guenther 
(1994). We can use this spectroscopy, together with photom- 
etry which exists for this star, to work through the steps we 
outline above and thus determine the appropriate parame- 
ters. It should be noted that the spectroscopy and photome- 
try are not simultaneous, and further observations would be 
necessary to fully vindicate the model. Guenther's spectra 
were not flux calibrated, so the flux variability needed to 
calculate Z is not directly available. Line fluxes should ide- 
ally be used for this measurement, however, Guenther notes 
that the variations in the equivalent widths of these lines for 
DG Tau are mostly due to variations in the line strengths, 
rather than continuum variations. In view of these consid- 
erations, we feel justified in using equivalent widths here to 
calculate Z. The approximate peak to peak variations of the 
lines studied by Guenther are shown in Table 1 

We will neglect Ca II 8662 and H a as they are lower 
energy, and less likely to be purely accretion driven. Taking 
a simple average of the other lines gives us Z = 1.96. From 
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20 


B 
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1.45 


55 


23 


70 


V 


0.35 


1.38 


62 


33 


91 


R 


0.3 


1.32 


68 


79 


100 


I 


0.22 


1.22 


78 







Table 2. Calculation of underlying photospheric and disc fluxes 
in each band for the T Tauri star DG Tau, and comparison with 
values estimated from an SED published by Basri & Bertout 
(1989) and photospheric templates for M0 and K7 dwarfs. Col- 
umn two shows the amplitude of the variability in each band, col- 
umn three the corresponding ratio of maximum to minimum flux. 
In column four we show the fraction of the flux attributable to 
underlying photosphere plus disk flux in Basri & Bertout 's SED. 
Columns five and six show the values obtained using our calcula- 
tions, column five assuming an M0 photosphere, and column six 
a K7 photosphere. 

equation 3 assuming 7 = 2, we can deduce that this implies 
an angle of 34 degrees between the axis of the dipole field 
and the plane of the disc. Photometry of DG Tauri taken 
by Bouvier et al (1993) shows periodicity in U,B,V,R and 
I bands. Applying the equation above to the amplitudes of 
the periodic variations found by them and using Z = 1.96 
gives the results shown in Table 2 Also shown are the values 
obtained by comparing a spectrum of DG Tau from a paper 
by Basri & Bertout (1989) to an M0 photosphere, as used 
by them, and to a K7 photosphere, also taken from their 
paper. 

It was not possible to obtain figures for the I band, as 
the published data do not extend far enough into the red. It 
can be seen from Table 2 that our values lie between those 
for the two photospheres. We therefore conclude that our 
method predicts that the photosphere of DG Tau is some- 
where between an M0 and a K7, possibly a K8. 



5 CONCLUSIONS 

We have discussed some consequences of the recently pro- 
posed beat frequency modulation model of T Tauri accre- 
tion, and showed how these consequences can be used to 
determine certain parameters of individual systems. 

The method outlined above provides a measurement of 
the underlying stellar and disc fluxes, independent of any 
spectral energy distribution (SED) modelling. This in itself 
is a very useful measurement, as it allows more tests to be 
applied to such SED calculations, and will also help to con- 
strain the positions of T Tauri stars on the HR diagram. This 
is perhaps most important in determining T Tauri ages. 

The model in this paper presents the first measurement 
of the relative orientation of the axis of the magnetic field 
and the disk in a T Tauri system. As illustrated in Section 
3, this orientation may have a profound effect on the inner 
structure of the accretion disk and infrared variability. These 
orientation diven effects, and their observables are currently 
a topic of further research and will be presented at a later 
date. 

Magnetically controlled accretion in T Tauri Stars has 
recently received much supporting evidence, both observa- 
tionally and as a model to explain their slow rotation rates 
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(Edwards 1995). It has also been hypothesized that the stel- 
lar magnetic field may be relevant in driving a stellar wind 
or jet (Edwards 1995; Shu et al 1994). In all of these models, 
the orientation of the magnetic field is assumed, for simplic- 
ity, to be perpendicular to the surrounding disc. The positive 
identification of a T Tauri Star undergoing periodic accre- 
tion at the beat frequency of the star and disc will yield the 
orientation of the magnetic field relative to the surrounding 
disc. This will be important not only as a constraint on the 
models of magnetospheric accretion and their possible rela- 
tion to stellar winds but also as a constraint on the physics 
of magnetic fields in pre-main sequence stars. 

The search for periodicities in emission line strengths 
is still in its early stages. These lines are strongly linked 
with other accretion signatures, however, and so it is to be 
expected that they will have such periodicities in common 
with the continuum. Recent studies appear to bear out thois 
expectation (Guenther 1994) A spectroscopic campaign is 
currently being planned which will as part of its aim, at- 
tempt a measurement of Z for a small number of T Tauri 
stars. 
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